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Abstract 

The Anderson model of local-state conduction electron mixing is applied 
to the problem of interlayer magnetic coupling in metallic multilayered 
structures with palladium (Pd) spacer layers. An oscillation period of 5 
spacer monolayers and the tendency towards ferromagnetic bias of the 
interlayer magnetic coupling that we obtain are consistent with the ex- 



Q\ ! perimental data. 

wo: 

The discovery of oscillatating interlayer magnetic couplings between ferromagnetic layers 
separated by a nonmagnetic metallic spacer JTJ and of the related giant magnetoresistance 
effect 0, has stimulated a lot of experimental and theoretical activity. It has been shown 
that the periods of the coupling are related to the topology of the Fermi surface of the spacer 
layers. This interpretation has been confirmed by model and first-principles calculations, 
and is also supported by experiments || |J. There are, however, other aspects of the 
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coupling, e. g., the bias (ferro- or antiferro- magnetic) of the interlayer magnetic coupling, 



. which have not been fully explained. 

For Fe(OOl) layers separated by Pd(OOl) spacers of thickness between 4 and 12 ML 
the interlayer magnetic coupling is observed to have a strong ferromagnetic bias as seen 
in the experiments ||. Above a 13 ML thickness of the Pd spacer the coupling begins 
to be antiferromagnetic. Metallic Pd is believed to be near the threshold of becoming 
ferromagnetic. The non-relativistic calculations of Moruzzi and Marcus || and of Chen et 
al. [[7] predicted the onset of ferromagnetism in fee palladium with a 5% expanded lattice. In 
a recent publication the ferromagnetic bias of the coupling in magnetic multilayer structures 
with a Pd spacer is explained in terms of the Pd as an almost ferromagnetic media [Ffl. 
Alternatively, in this paper, we interpret this ferromagnetic bias to be a consequence of a 
competition between RKKY-like and superexchange couplings, with RKKY coupling being 
dominant. 

The RKKY-like coupling comes from intermediate states which correspond to spin exci- 
tations of the Fermi sea. States corresponding to electron-hole pair production in the Fermi 
sea, with an attendant spin-flip, contribute to the RKKY coupling as ||; 

jRKKY(q)= 2J [ —r —y 2 — + «.] , (1) 

ni,n2, fc \ t n2k' t + ) fcn 2 fc' t nik 

where 6 is a step function, e F is the Fermi energy, k' = k + q + G, G is a vector of 
the reciprocal lattice, e + is the energy of the local impurity state, and V n k represents the 
strength of the s — d mixing interaction |J . 
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Figure 1: RKKY-like coupling, Jrkky{(Iz), and superexchange coupling, js(c[z), for 
Pd(OOl) spacers in reciprocal space calculated with E h = 0.08 Ry. V is expressed in 
rydbergs. 



The superexchange coupling arises from charge excitations in which electrons from local 
states are promoted above the Fermi sea (one from each layer) providing a second contri- 
bution to the coupling in addition to the RKKY coupling ||: 

/ \ r I Vnik | 2 I Ki 2 fc' | 2 ^(^nife ~ ^F)0(e n . 2 k' ~ Sp) 

Mq) = - ;L i—fZ —~Y2 — + c.c.]. (2) 

ni,n 2 ,k \ t n 2 k' t + j t-nifc fc + 

The real space coupling between two sheets of spins can be obtained by Fourier trans- 
forming Eqs. (1) and (2) ||, with the coupling in the real space given by, 

J t{ z ) = 7T / dq z j(q z ) cos(q z z) , (3) 

Z7T JO 

where a is a lattice constant, and z is in the direction perpendicular to the magnetic layers. 
The sign is chosen so that positive Jt(z) signifies ferromagnetic coupling. 

Using the Slater-Koster parameters [IIJ , one can easily diagonalize small matrices (9x9 
for a typical transition metal) to obtain the energy bands and density of states (DOS) for 
fee Pd. The electron wave function | n, k > is a Bloch state belonging to band n and wave 
vector k, and is expressed as linear combinations of localized orbitals: 

I n, k >= J= £ e ik - R " ]T a m (k) Ui (v - R u ) , (4) 
viV „ 



where iV is the number of cells in the material considered, H u is a lattice vector, m(r — H u ) 
is the ith orbital basis function, and «m(k) is a (real) normalized eigenvector compo- 
nent determined by diagonalization of the single-particle Hamiltonian. We use a plausi- 
ble approximation V ni kV* 2k , = V 2 M* ik n2k , (q) ||, where the matrix element is defined as 
M ni k,n 2 k' — < n ik I e * qr I n 2^ >■ The explicit expression for the matrix element is 



M ni k,n 2 k' — J! 



e lkn " E a ni (k) a n2j (k') J dv u t (r) e iq ' r Uj (r - R u ) . (5) 



The essential conditions for the simplication of this matrix was already discussed by Call- 
away et al. [pTlfl . Mj(r) are approximated as Clementi wave functions for the d states 
and a„j(k) can be related to the Slater-Koster parameters in Ref. [10]. 

We consider one local level below Ef for simplicity and set e + = Ef — Eh, where Eh 
is the energy required to promote an electron from an occupied local magnetic impurity 
level to the Fermi level. Based on the band structure of bulk paramagnetic Pd, we have 
calculated the couplings 3rkky{<Iz) and jsilz) with E h = 0.08 Ry, as shown in Fig. 1. 

The couplings in real space are plottted in Fig. 2. The dashed line, dotted line and 
solid line are for RKKY-like, superexchange, and RKKY + Superexchange couplings, re- 
spectively. We see that the superexchange interaction gives a small contribution to the 
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Figure 2: Interlayer magnetic coupling in the real space. The dashed line and dotted line 
are for RKKY-like and superexchange couplings, respectively. The solid line is for the total 
coupling (RKKY+superexchange). 



Figure 3: Density of states for Pd. 



coupling, and the total coupling has a strong ferromagnetic bias. This tendency for a 
ferromagnetic bias resembles the experimental observation in Fe/Pd(001) trilayered struc- 
tures 0[. The 5 ML oscillation period in the calculated interlayer magnetic coupling J(z), 
as shown in Fig. 2, corresponds to the peak at q z ~ 0.4 — in Jbkky (Qz)- It agrees with the 
experimental period of 4 — 5 ML in Fe/Pd/Fe(001) for trilayered structures ||. 

We explain the result of ferromagnetic bias for multilayers with Pd spacers as being 
due to the structure of the Pd DOS and the location of the Fermi level, as shown in 
Fig. 3. In particular, the fact that the Fermi level falls above a peak in the DOS followed 
by a relatively smooth and structureless DOS, results in a relatively small superexchange 
contribution above ~ 4 ML, and leads to a ferromagnetic bias driven by the large RKKY 
coupling. 

To confirm our explanation of the cause of the ferromagnetic bias we use a free-electron 
gas model which enables us to obtain analytic results for the couplings. In a previous 
study [[T2J , we showed that in the free-electron gas approximation, RKKY + Superexchange 
coupling resembles pure RKKY coupling, but without any magnetic bias. To illustrate the 
effect of a peak in the density of states on top of a free electron-like background, we use a 
"toy model" calculation by adding a Lorentzian shaped peak to the DOS of the free electron 
gas, 

where the position of the peak is at E p = p ep, and h adjusts height of the peak (small h 
corresponds to a large peak). For example, by fixing the position (p = 0.9, below the Fermi 
level) and increasing the height of the peak (e.g., h = 0.3), ferromagnetic bias occurs in 
the coupling, as shown in Fig. 4. We noted that in multilayer structures with a Cr spacer, 
RKKY + superexchange coupling gives an antiferromagnetic bias due to the structure of 



the DOS, with a peak above, but near to the Fermi level ep [12]. This is also confirmed in 
our "toy model" calculations. 

The coupling observed in the experiments contains a bilinear exchange interaction J\ 
and a biquadratic exchange interaction J 2 . In conventional notation, J exp — Jy — 2J 2 [El. 
Here, RKKY-like and superexchange couplings are contained in the bilinear coupling, J\, 



Figure 4: The total interlayer magnetic coupling for a densities of states with a peak 
below, but near, the Fermi surface. Ferromagnetic bias appears in the coupling. 
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Figure 5: Interlayer magnetic coupling as a function of Pd spacer thickness. The solid line 
and filled circles are our calculations and the experim ental data observed at T = 77 K Jl3| , 
respectively. The theoretical (experimental) results are refered to the left (right) scales. 



and a positive J 2 favors a perpendicular magnetic coupling. The ratio of J%/ J\ observed 
in the multilayer structures with Cr spacers is about 0.3 - 0.5 [[R|. Magnetic multilayer 



structures with Pd spacers also have a relatively large biquadratic exchange interaction 
J 2 @- With a proper choice of positive J 2 (J2/ Ji ~ 0.5 at 13 ML), the bias can switch to 
antiferromagnetic for spacer thickness greater than 13 ML, as is observed experimentally. 
One can see that our calculated interlayer magnetic coupling with Pd spacers can be used 
to explain the experimental data as shown in Fig. 5. 

In summary, our model calculation has been able to reproduce the two salient features 
of the interlayer magnetic coupling in Fe/Pd(001) multilayer structures: large but rapidly 
decreasing ferromagnetic bias, and a 5 ML oscillation period. The ferromagnetic bias 
arises from the competition between the RKKY-like and superexchange couplings due to 
the special features of the palladium DOS: relatively large peak below, but near to the 
Fermi level, and a small DOS above the Fermi level. 
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